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We show that the entanglement witnesses based on local orthogonal observables which are intro- 
duced in [S. Yu and N.-L. Liu, Phys. Rev. Lett. 95, 150504 (2005)] and [O. Giihne, M. Mechler, 
G. Toth and P. Adam, Phys. Rev. A 74, 010301 (R) (2006)] in linear and nonlinear forms can be 
optimized, respectively. As applications, we calculate the optimal nonlinear witnesses of pure bi- 
partite states and show a lower bound on the I-concurrence of bipartite higher dimensional systems 
with our method. 
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I. INTRODUCTION 

Entanglement is one of the most fascinating features of 
quantum mechanics, which has recently been recognized 
as a basic resource in quantum information processing 
such as teleportation, dense coding and quantum key dis- 
tribution [if, Thus, it becomes particularly important 
to detect and quantify entanglement [sj . Despite a great 
deal of effort in the past years, lots of things are still 
unclear to us in this field (see the reviews [J d, [1] and 
references therein). Nevertheless, on the one hand, sev- 
eral sufficient conditions for detection of entanglement 
have been found, such as the famous Peres-Horodecki 
positive partial transpose (PPT) criterion 0, [Ij , realign- 
ment criterion [9|, entanglement witnesses (EWs) 
local uncertainty relations (LURs) [ll|, [l^j. Bell type 
inequalities [H, O, [iBj, etc. PPT criterion is necessary 
and sufficient for 2 x 2 and 2x3 systems, but only neces- 
sary for higher dimensional cases Q. It is believed that 
realignment criterion complements PPT criterion since 
it can detect many entangled states which PPT crite- 
rion cannot detect. More easier way to detect entangle- 
ment experimentally is using EWs^which have recently 
been generalized to nonlinear EWs \vf\ . On the other 
hand, a considerable amount of effort on quantification 
of entanglement has also been made. For instance, Woot- 
ters has analytically derived a perfect measure of 2 qubits 
[l8| . which is so-called concurrence. Furthermore, gener- 
alized concurrence in bipartite higher dimensional cases 
[Tgl . [20| . such as I-concurrence f2o!], has been pointed 
out as well. Unfortunately, the I-concurrence of mixed 
states is given as a convex roof for all possible ensem- 
ble realization. Therefore, it is generally difficult to be 
calculated. Lately, lower bounds on I-concurrcnce have 
attracted much interest [2l|, [12, [H, [13] , which are rela- 
tively easier than I-concurrence itself to get. 

Recently, Yu and Liu have introduced an entangle- 
ment witness [Eq. ([3])] based on local orthogonal ob- 
servables (LOOs) in Ref. [251]. Moreover, Giihne et 



al. have generalized the witness to the nonlinear form 
[Eq. (HI] via local uncertainty relations 26]. Both of 
the witnesses have a common property that each set of 
LOOs in the witnesses can be replaced by any other 
complete set of LOOs, thus one does not know which 
set of LOOs is the best one for the witnesses. Actu- 
ally, the witnesses using different set of LOOs can ob- 
tain distinct results. For example, the Bell state (JOO) -I- 
\ll))/y/2 can be detected as entangled states by the 
liner witness under the set of LOOs: {ctj,, o-y , CTz, / j'^/v^, 
{(Ta:, — (Tj,, (T^ , /}^/-\/2, but cannot be detected under the 
LOOs: {a,,ay,a„I}'^/V2, {cj^,ay,a,,I}^ 1^2. There- 
fore, it is necessary to investigate the optimal case. In 
this paper, the optimal witnesses for the linear and non- 
linear forms will be presented. As applications, we will 
calculate the optimal witnesses of pure bipartite states 
and show a lower bound on the I-concurrence of bipar- 
tite higher dimensional systems. 

The paper is organized as follows: Sec. II presents the 
optimal witnesses of linear and nonlinear forms, which 
are constructed by LOOs. In Sec. Ill we calculate the 
optimal nonlinear witnesses of pure bipartite states based 
on our method. Moreover, we obtain a lower bound of I- 
concurrence in bipartite systems. Sec. IV discusses what 
happens if the dimensions of the subsystems A and B are 
not the same. 



II. OPTIMAL WITNESSES BASED ON LOOS 

For convenience, we consider a, d x d bipartite system, 
just as Refs. [H, [2^ did ( in Sec. IV we will discuss 
the situation when dimensions of subsystems A and B 
are not the same). Each subsystem has a complete set 
of local orthogonal bases {G^} and {Gf }, which are so- 
Such a basis consists of d^ observables and 



called LOOs. 
satisfies: 



MGtGf) = Tr(Gf Gf ) 



6ki 



(1) 



Any other complete set of LOOs relate to the original 
one by an orthogonal d^ x real matrix, i.e., 
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where OO'^ = 0^0 = O'O'^ = O'^C = /. 

In Ref. [2^, a linear witness was introduced as fol- 
lows (for convenience, the witness has been written in an 
equivalent form introduced in [26j ). 



(3) 



where {G^} and {Gf } are arbitrary complete sets of 
LOOs for subsystems A and B. Later, Ref. [2^ provided 
a nonlinear form. 



^Y^iGt^i-mG^y. (4) 



For every separable state p, it must satisfy that TrWp > 
and J-{p) > 0. Conversely, if any state violates one of 
the two inequalities, it is entangled indeed. 

In Refs. [2^, [2^, there is a little mention involving 
that how to choose a set of LOOs so that TrWp or J^{p) 
gets its minimum, and obviously the minimum means a 
optimal one, since one can obtain distinct results by us- 
ing different sets of LOOs. Consider the simple example 
\^P+) = (|00) + |ll))/\/2 introduced in Sec. I. Under the 
set of LOOs {ax,cry,az,I}^/V2, {ax,(Jy,(JzJ}^/V^, 
Tr(W|i/'+)(V'+|) = and = 0, with 

which one cannot conclude that is entangled. 

However, under the set of LOOs {da,, (Ty, cr^, /}"^/\/2, 
Tr(W|^+)(^+|) = -1 and 
T{\'4>^){ip^\) = —1. It suggests that |?/^+) has entangle- 
ment. Therefore, it is meaningful to obtain the minimal 
one. In the following, we will show that the minimum is 
invariant under local unitary (LU) transformations, and 
obtain an analytical formula of the minimum. 

Lemma 1. For a given state p, the minimum of TrWp 
[J-{p)] is LU invariant. 

Proof — (Reductio ad absurdum) For a given state p, 
suppose that under the set of LOOs {M^}, {Mf } TrW/o 
gets its minimum Li. We operate an arbitrary LU 
transformation to p, i.e., p' — Ua® UbpU\ ® Ug. For 
the state p' , suppose that under the set of LOOs {M^}, 

{Mf } TrWp [^(p)] gets its minimum L2. 

Case i. Li > L2. For the state p, under the set of 

LOOs {u\m^Ua}, {uIm^Ub}, TrWp [J'ip)] is equal 
to L2. It is a contradiction to that Li is the minimum of 
TrWp [J'ip)]. 

Case ii. Li < L2. For the state p' , under the set of 
LOOs {UaM^u\}, {UbM^uI}, TrWp' [J^(p')] is equal 
to Li. It is a contradiction to that L2 is the minimum of 
TrWp' [T{p% 

In a word, if Li ^ L2, a contradiction is derived im- 
mediately. Therefore, Li = L2 always holds and the 
minimum of TrWp [.?^(p)] is LU invariant. □ 

Remark.— From an experimental point of view, it is 
valuable for the minimum to satisfy LU invariant condi- 
tion, since a shared spatial reference frame is no longer 
needed when one makes a measure of the minimum 12711 . 



Theorem 1. The minimum of TrWp is equal to 1 — 
Sfc '^kiP'), where iTfc(/i) stands for the fcth singular value 
of real matrix p, which is defined as pi„i — Ti{pGf(^G^). 

Proof.— Before embarking on our proof, it is worth 
noticing that a similar result of Theorem 1 has also been 
pointed out in (25| . However, for a convenience to un- 
derstand Theorem 2, we insist on providing a complete 
proof. For a given state p, we choose an arbitrary com- 
plete set of LOOs {G^}, {Gf }. Define that 



Pi,„ = Tr(pGf ®G^), 
and the density matrix can be written as: 
P = Ep;mG/^(»G^. 

l,m 



(5) 



(6) 



According to Eq. any other complete set of LOOs 



{Gt}, {G|} can be written as G^ = J^i UkiG, 
J2m^kmG^, where U and V are (P x (P real orthogo- 



nal matrices. 
Therefore, 



i.e. UU' ^ U'U ^VV ^ V'V = I. 



minTr(Wp) = 1 - maxE(G^ «) Gf ) 

k 

= 1 - max J2 ^kiVkm {Gf ® G^) 

k Im 

= 1 - max E E Ukl VkmPhn 

k Im 

= 1 - maxE[C/^y^]fcfc 

k 

= 1 -maxTr(C/MV"^). (7) 

Moreover, 

maxTr([//iy^) = maxTr(/iy^[/) = ^ crfe(p), (8) 



where we have used the following theorem [28j: 

Let A e Mn he a given matrix, and let A = V^W^ be 
a singular value decomposition of A. Then the problem 
max{Re trAU : U € M„ is unitary} has the solution 
U = WV\ and the value of the maximum is cri{A) + 
• • • + (^n{A), where {<Ti{A)} is the set of singular values 
of A. 

Notice that p is a real matrix and its singular 
value decomposition can be written as p = U'^TiV, 
where U, V are real orthogonal matrices and S = 
diag{ai{p),a2{p), ■ ■ ■ ,o'^2(p)}. When U — U and V — 
V, Tr(J7pl/"^) gets its maximum 'Ylik'^kip)- In other 
words, under the new complete set of LOOs {Q^}^ 
where = Y^MuGt = E^VkmGg, W - 1 - 
Y.k ^k ®^k ^ TrWp gets its minimum 1 - Y,k '^k{p)- □ 

Remark.— In fact, it is equivalent to the realignment 
criterion when TrWp gets its minimum (25j . Note that 
under the new complete set of LOOs {^/^}, {^/^}, the 
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density matrix can be written in its operator-Schmidt 
decomposition form (29j : 



LOOs transformations, i.e. 



(9) 



The realignment criterion states that if p is separable 
the sum of all akifJ-) is smaller than 1. It is equivalent to 
minTrWp > 0. Hence, it is concluded that any entangled 
state detected by a witness of Eq. ([3]) must violate the 
realignment criterion. 

Example.— Let us consider a noisy singlet state intro- 
duced in Ref. [Ill, p = p|V's)(V's| + (1 - p)psep, where 
lips) stands for the singlet state (|01) — |10))/-\/2 and 
the separable noise is p^ep = 2/3|00)(00| + 1/3|01)(01|. 
Actually, the state is entangled for any p > [2^]. Un- 
der the complete set of LOOs {— c^, —cry, —crz,I}^/V^, 
{ax,(Ty,crz, I}^ /V^, the witness of Eq. ([3|) can detect 
the entanglement for all p > 0.4. However, the optimal 
witness using Theorem 1 can detect the entanglement 
for all p > 0.292, which is equivalent to the realignment 
criterion. 

Theorem 2. The minimum of J-{p) is equal to 1 — 
J2k (^kir) — (Trp^ -I- Trp^)/2, where crfe(r) stands for the 
kth singular value of matrix r defined as Tim = {Gf (8) 
Gg)-{Gf®I){I®Gf^). 

Proof.— For a given state p, we choose an arbitrary 
complete sets of LOOs {G^},{Gf}, and calculate the 
real matrix r according to the definition: 



Tim - {Gf <E> GZ) - {Gf ® /) (/ ® Gfj . (10) 



Similarly to Theorem 1, any other complete set of LOOs 

{G^}, {GQ can be written as = UuGf, Gf = 
VkmGm, where U and V are tP x (P real orthogo- 



nal matrices 
Therefore, 



i.e. UU' ^ U'U = VV' =.V'V = I. 



E 

k 



k W 



W 



U' U]u'{Gt®I){G^®I) 



I 



{Gt®iY 



= Trpi, 

where pA is the reduced density matrix after tracing over 
subsystem B. Without loss of generality, substituting Eq. 
([TT|) into ^Ylii^G^ (8i /)^, one can obtain the final result 

Trp2 . Similarly, Y^k^l ® G^f ^Y.M® Gff = Trp| 
holds. 

-^(|m)(n| + |7i)(m|) 1 < m < n < d, 
G^ = { -j=ii\m){n\-i\7i,){m\) 1 < m < n < d, (11) 



\/2 

|m) (m| 



1 < m < d. 



A\T 



Gf = (Gt) 



(12) 



where {Im)^} and {|to)b} are the standard complete 
bases. Thus, 

max[5](G^ ® Gf ) + i 5](G^ ® / - / ® Gf )2] 

k k 

= lj2^{Gt®lf + {I®GSf] 

k 

+max[5]((G^ ® Gf ) - (G^ ® /) (/ ® Gf ))] 
fe 

= ^ (Trpi + Trp| ) + max ^ ^ Uki Vkm Tim 

k hn 

= i (Trp^ + Trpl ) + max ^ [C/rV ^] 



i(Trpi + Trp|)+^afe(r). 



(13) 



min[l - J2{G^ ^G^)--Y,{G^^I-I^G 



= 1 max[^(G^ ® Gf ) + - ^(G^ ® / - / ® G^)']. 

k k 

Moreover, 

Y^iGt^I-I^GSy 

k 

= J2i(G^^^I)' + {I^Gff-2{Gl^I){I^Gf)l 

k 

where J^kiG^'^I)'^ and J^ki-^^G^)^ are invariant under 



In other words, minJF(/9) = 1 — J2k'^k{T) — (Trp^ -|- 
Trp|)/2. ^ □ 

Example.— Bennett et al. introduced a 3 x 3 bound 
entangled state constructed from unextendible product 
bases in Ref. [soj : 

iV'o) = ^I0>(|0) - |1)), IV^i) = ^(|0) - |1))|2), 

\^,) = ^|2)(|1) - |2)), - - |2))|0), 

|V'4) = ^(|0) + |1) + |2))(|0) + |1) + |2)), 



1 4 



(14) 
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Let us consider a mixture of this state with white noise, 

p{p)=PP+{l-p)^. (15) 

Using the reahgnment criterion, one finds that the state 
p(p ) stiU has entanglement when p > 0.8897. In Ref. 
[26|, it is found that the state p{p) must be entangled 
for p > piur — 0.8885 using the nonlinear witness Eq. 
(HI) (but not the optimal one). According to Theorem 2, 
one can obtain an optimal witness of Eq. ([4]) and find 
that when p > popt = 0.8822 the state is still entangled. 
Obviously, the optimal witness is stronger than the one in 
Ref. [2^. In addition, in Sec. Ill we will present a lower 
bound on I-concurrence for the state based on Theorem 
2 (see Fig. 1). From the figure, it is worth noticing that 
the bound is positive when p > popt — 0.8822. 



Therefore, 

(18) 

k m<n m<n 

i(Trp^ + Trp|)-^/i2, (19) 

i 

CmaAm = iY.^^^- (20) 

i 

Note that Eq. ([^0]) has also been derived with another 
totally different method in Ref. [23*1, and it completely 
accords with our result. Compared with the method in 
Ref. [2^, Theorem 2 in this paper is more general, i.e., it 
suits not only bipartite pure states but also any bipartite 
mixed state. 



B. Lower bound on the I-concurrence 



III. APPLICATIONS 

In this section, the optimal nonlinear witnesses of pure 
bipartite states will be obtained using Theorem 2. More- 
over, we will show a lower bound on the I-concurrence 
of bipartite systems by means of our method. Be- 
fore embarking on our investigation, we first define that 



)Gf), and ob- 



viously Cmax = J2k '^kir) + i'Trp\ + Tr/7|)/2 according 
to Theorem 2. 



A. Optimal witnesses of bipartite pure states 

Let us calculate Cmax of a bipartite pure state {ip) with 
its Schmidt decomposition \ip) = J2i v^N)- 

Since Schmidt decomposition of a pure state is a LU 
transformation, Cmaxi\'4')) remains invariant after the 
transformation according to Lemma 1. Therefore, we can 
directly use the Schmidt decomposition form for conve- 
nience. We choose a complete set of LOOs Eq. ([Tl]) and 
Eq. ((T2)) for A and B subsystems, respectively (obviously 
any other complete set of LOOs can be chosen and it does 
not affect the final result). 

According to Theorem 2, 

Tim = {Gt®G^)~{Gt®I){I®G^:, 

= [D®D®T]im, (16) 

where D = diag{^pip2, ■ ■ ■ , ^/PmPn • ■ ■ , y/Pd~iPd} and 



T = 



f Hi -pi ■■■ -P'lPd \ 

~PlP-2 - ■ ■ • -P'2Pd 



\ -piPd -P2Pd ■ ■ ■ Pd- pI J 



(17) 



I-concurrence of a bipartite pure state is given by 
C'dV')) = ■\/2(l — Trp^), where the reduced density ma- 
trix Pa is obtained by tracing over the subsystem B. It 
can be extended to mixed states p by the convex roof, 

for all possible decomposition into pure states, where 
Pi > and J2zPi = 1- 

Several bounds have already been derived [21], [H, [H, 
[2^ . e.g., an analytical lower bound based on PPT cri- 
terion and the realignment criterion has been shown by 
Chen et al, 



C{p)> 



m{m — 1) 



(max(||p^-||,|l7^(p)||)-l), (22) 



where Ta, TZ and || • || stand for partial transpose, re- 
alignment and the trace norm (i.e. the sum of the sin- 
gular values), respectively. In Ref. [11], another bound 
based on LOOs has been obtained, which has used Eq. 
(HOI) and the fact that Y.iPi^rnaxHi'i)) > J^iPi^ili'i)) > 
C{J2iPi\'^i){^i\), (for convenience, the lower bound has 
been rewritten in an equivalent form) 



C{p)> 



m{m — 1) 



(£-1). 



(23) 



Notice that Eq. (|23p holds for arbitrary set of LOOs, 
including the optimal one. Therefore, a tighter form of 
Eq. can be obtained according to Theorem 2, 



Cip)> 



m{m — 1) 



1), 



(24) 



where Cmax = T,k ^k ir) + (Trp\ + Trp|)/2. Since the 
entanglement criteria based on LURs are strictly stronger 
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Theorem 1 and Theorem 2, respectively. 
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FIG. 1: Two lower bounds of I-concurrence for the state p(p). 
One is the lower bound based on realignment criterion (dashed 
line), the other is obtained from jCmax (sohd line). 



than the realignment criterion [26| , the following inequal- 
ity can be concluded. 



C{p)> 



m{m — 1) 



{max{\\p'^^\\,C,naxip))-l). (25) 



For example, reconsider the bound entangled state Eq. 
(fT4|) . Because it belongs to PPT entangled state, the 
lower bound based on PPT criterion is unhelpful. One 
can obtain that C{p) > 0.050 via the realignment crite- 
rion, and C{p) > 0.052 has been gotten in Ref. [11] by 
using Eq. (1^ . In fact, Cmax{p) can be directly calcu- 
lated, and it suggests that C{p) > 0.055 via Eq. ([M]) . 
which is better than the one in Ref. [l^. Furthermore, 
one can consider the bound entangled state with white 
noise, i.e. Eq. ([TS]). The lower bounds of I-concurrence 
for p{p) have been shown in Fig. 1. Therefore, the lower 
bound based on Cmax has been strictly improved com- 
pared with the one based on the realignment criterion 
and provided a tighter form of Eq. 



IV. DISCUSSION AND CONCLUSION 

During the last two sections, we consider a simple sit- 
uation: the 0? X c? bipartite system for convenience. How- 
ever, if the dimensions of the Hilbert spaces Ha and Hb 
are not the same, what will happen? Actually, it does 
not affect any one of the conclusions in Sec. II and Sec. 
III. 

Without loss of generality, suppose that m — 
dim(7i/i) < n ~ dim(7YB). There are m? elements in 
a complete set of LOOs {G^}, and in? elements in {G^}. 
Therefore, we need to reconsider Eq. ([5]) and Eq. (fT5|) in 



max^ ^ UkiVumPim = maxTr(C/^y'^), (26) 

k Ira 

max^ ^ UuVkmTinr = maxTr(f/rl^^), (27) 



k In 



where U is an 
T are x 



■? real orthogonal matrix; p, and 
real matrices; V belongs to v? x real 
orthogonal matrices. The two equations have the same 
form, so we .just need to consider Eq. (P?]) for instance. 
As Ref. \2A did, one can define that = for A; = 



TO 
^2 



^ did, one can define that 

-|- 1, • • • , n^. Thus, the matrix r is changed into an 
X T? real matrix, i.e.. 



(28) 



where stands for an (r? — rr?) x r? matrix with every 
element being equal to 0. 

Define that U' = U®I, where / is an {n^ — m?) x [n^ — 
TO^) identity matrix. It is easy to see that U' is an n? x 
real orthogonal matrix since U belongs to m? x m? 
orthogonal matrices. 

Notice that {I = 'n? — m^) 



//x; 



0;xn2 



0/xn2 



which means that Tr[t/'r'V"^] = TrprV'^]. Therefore, 

maxTr[C/Ty'^] = maxTr[C/'T'y'^] 
= maxTr[r'y^C/'] 

^a,(r'). (29) 



Since t't''^ = 1''"'""^] ®Oixi, t't''^ and tt'^ have the same 
nonzero eigenvalues. Hence, 



^a,(r')-^a,(r). 



(30) 



Consequently, Eq. ^ and Eq. ^ suggest that The- 
orem 1 and Theorem 2 still hold even if the dimensions 
of subsystems A and B are not the same, and the appli- 
cations in Sec. HI which have used the Theorem 2 can 
also be extended to this case. 

In conclusion, we have optimized the linear and the 
nonlinear entanglement witnesses based on local orthog- 
onal observables, which are introduced by Yu, Liu and 
Giihne et al. respectively, and several examples have 
been given as well. Moreover, we have obtained the opti- 
mal witnesses based on LOOs in pure bipartite systems 
and a lower bound on the I-concurrence of bipartite sys- 
tems as applications of our method. In fact, Theorem 2 
presents a separability criterion with Ky Fan norm of r, 
the covariance term defined in p7j . Similarly, another 
separability criterion with Ky Fan norm of correlation 
matrix has been shown in [3l|. It is worth investigating 
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deeper relation between this two criterions. In addition, 
the 'optimal' in this paper is in the sense of choosing the 
best complete set of LOOs such that the witness gets 
its minimum, which has little relation with traditional 
optimal EWs [32,]. 

Note added. Recently a similar result has been shown 
in (33] , which is based on covariance matrix criterion. 
Interestingly, Proposition 3 in can be optimized to a 
similar form as Theorem 2 in this paper. 
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